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1 Field Of The Invention 

2 The invention relates generally to the field of computer graphics, and more specifically to 

3 generation of meshes representing surfaces of objects. The invention specifically provides a 

4 computer graphics system and computer-implemented method for generating a coarse-level mesh 

5 from a fine-level mesh, which may be used in multi-resolution surface editing. 

6 Background Of The Invention 

7 Generally, in computer graphics, objects are represented as surfaces, with the surfaces being 
represented by meshes. A mesh consists of a set of vertices, or points in three-dimensional space, 

—9 interconnected by edges. The edges define polygonal faces, which may be in the form of triangles, 

ill 

rHO quadrilaterals, and so forth. In some computer graphic operations, it is desirable to generate a 

f=f 1 representation of a surface at a finer resolution than a current representation. There are several 

U I 2 popular methodologies for generating a representation of a surface at a finer resolution than a current 

5 

□ 3 representation, including a Catmull-Clark surface subdivision methodology and a Loop surface 

f fi 

p| 4 subdivision methodology. In other operations, it can also be desirable to generate from a mesh 

^1 5 representing a surface at one level of resolution, a mesh at a coarser resolution. A methodology for 

Ll 6 generating a representation of a surface at a coarser resolution than a current representation, and a 

17 methodology for generating such a surface representation, based on the Loop methodology, is 

18 described in D. Zorin, et al., "Interactive Multiresolution Mesh Editing, SIGGRAPH 1997, 

19 Computer Graphics Proceedings, pp. 259-268, particularly p. 262. 

20 The methodology described in Zorin relies on a method for smoothing triangular meshes that 

2 1 is described in G. Taubin, "A signal processing approach to fair surface design," SIGGRAPH 1 995, 

22 Computer Graphics Conference Proceedings, pp. 351-358. Taubin's smoothing methodology 

23 consists of two linear filtering steps, the first being a low-pass filter step, followed by a high-pass 

24 filter step. The high-pass filter step is used to cancel out a mesh-shrinking effect that is induced by 

25 the low-pass filter step. The filtering operations described by Taubin have the problem that, as a 

26 consequence of the linearity of the filtering procedure, artifacts referred to as Gibbs ripples are 
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1 introduced in the smoothed mesh as a result of discontinuities in the input mesh. When used in 

2 generating a coarse level mesh, Taubin's filtering produces Gibbs ripples at the places on the surface 

3 where the fine level mesh has sharp variations. Although Taubin's smoothing methodology may be 

4 configured with the help of parameters, no filter can play the role of an inverse of the Loop 

5 subdivision rule. Accordingly, in a multi-resolution representation of a Loop subdivision surface, 

6 there will be redundant detail information. 

7 Summary Of The Invention 

8 The invention provides a new and improved system and method for generating a coarse-level 
f s% 9 mesh from a fine-level mesh, which may be used in multi-resolution surface editing. 

CfflO In brief summary, the invention provides an arrangement for generating a coarse ("j") level 

r|l mesh representing a surface, from a finer ("j+1") level mesh surface representation. The 

ss. 

JI2 arrangement includes an indicator value generator and a coarse level mesh generator. The indicator 

013 value generator, for respective ones of the points in the finer level mesh surface representation, 

t z \4 evaluates an indicator function, the value indicating whether a subdivision-inverse filter 

::| 5 methodology or a least-squares optimization methodology is to be used to determine a position for 

ul6 the corresponding point in the coarse level mesh representation. The coarse level mesh generator 

H7 determines, for each of the points that is to be provided in the coarse level mesh representation, a 

1 8 position in response to the position of the corresponding point in the finer level mesh representation, 

19 in accordance with the one of the subdivision-inverse filter methodology or least-squares 

20 optimization methodology indicated by the indicator value generated by the indicator value 

21 generator. 

22 Brief Description Of The Drawings 

23 This invention is pointed out with particularity in the appended claims. The above and 

24 further advantages of this invention may be better understood by referring to the following 

25 description taken in conjunction with the accompanying drawings, in which: 
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1 FIG. 1 depicts a computer graphics system including an arrangement for generating a coarse 

2 level mesh representation of a surface from a finer level mesh of a surface, constructed in accordance 

3 with the invention; 

4 FIG. 2 depicts a mesh representing a surface to which the Loop surface subdivision 

5 methodology has been applied; 

6 FIG. 3 depicts a mesh representing a surface to which the Catmull-Clark surface subdivision 

7 methodology has been applied; 

8 FIGS. 4 and 5 depict diagrams useful in understanding operations performed by the fine-to- 

9 coarse level mesh generating arrangement; 

y]|0 FIG. 6 depicts a flow chart describing operations performed by the fine-to-coarse level mesh 

I pfl generating arrangement in connection with triangular meshes; and 

fij 

FIG. 7 depicts a flow chart describing operations performed by the fme-to-coarse level mesh 

063 generating arrangement in connection with quadrilateral meshes. 

c:? 

W| 4 Detailed Description of an Illustrative Embodiment 

•S3? 

5 FIG. 1 depicts a computer graphics system 10 including an arrangement for generating a 

16 coarse level mesh representation of a surface from a finer level mesh of a surface, constructed in 

1 7 accordance with the invention. With reference to FIG. 1 , the computer graphics system includes a 

18 processor module 11, one or more operator input devices 12 and one or more display devices 13. 

1 9 The display device(s) 1 3 will typically comprise a frame buffer, video display terminal or the like, 

20 which will display information in textual and/or graphical form on a display screen to the operator. 

21 The operator input devices 12 for a computer graphics system 10 will typically include a pen 14 

22 which is typically used in conjunction with a digitizing tablet 15, and a trackball or mouse device 

23 16. Generally, the pen 14 and digitizing tablet will be used by the operator in several modes. In one 

24 mode, the pen 14 and digitizing tablet are used to provide updated shading information to the 

25 computer graphics system. In other modes, the pen and digitizing tablet are used by the operator to 

26 input conventional computer graphics information, such as line drawing for, for example, surface 
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1 trimming and other information, to the computer graphics system 10, thereby to enable the system 

2 1 0 to perform conventional computer graphics operations. The trackball or mouse device 1 6 can be 

3 used to move a cursor or pointer over the screen to particular points in the image at which the 

4 operator can provide input with the pen and digitizing tablet. The computer graphics system 1 0 may 

5 also include a keyboard (not shown) which the operator can use to provide textual input to the 

6 system 10. 

7 The processor module 1 1 generally includes a processor, which may be in the form of one 

8 or more microprocessors, a main memory, and will generally include one a mass storage subsystem 

9 including one or more disk storage devices. The memory and disk storage devices will generally 
1 0 store data and programs (collectively, , 'information ,, ) to be processed by the processor, and will store 
|=| processed data which has been generated by the processor. The processor module includes 

connections to the operator input device(s) 12 and the display device(s) 13, and will receive 

Ml information input by the operator through the operator input device(s) 12, process the input 

J?4 information, store the processed information in the memory and/or mass storage subsystem. In 

Ijf addition, the processor module can provide video display information, which can form part of the 

1 6 information obtained from the memory and disk storage device as well as processed data generated 

iff thereby, to the display device(s) for display to the operator. The processor module 1 1 may also 

S include connections (not shown) to hardcopy output devices such as printers for facilitating the 

ill 

© generation of hardcopy output, modems and/or network interfaces (also not shown) for connecting 

lb the system 10 to the public telephony system and/or in a computer network for facilitating the 

2 1 transfer of information, and the like. 

22 The invention provides an arrangement for generating a coarse-level mesh from a finer level 

23 mesh. In the following, the fine level mesh will continue to be identified by index "j+1," and the 

24 coarse level mesh will continued to be identified by index "j." Before describing the operations 

25 performed by the fine-to-coarse mesh generating arrangement in connection with generating a 

26 coarse level mesh from a current representation, it would be helpful to describe two surface 

27 subdivision methodologies, namely the aforementioned Loop surface subdivision methodology and 

28 a Catmull-Clark surface subdivision methodology, which are used to generate a finer ( M j+1) level 

29 mesh from a coarse ("j") level mesh. The invention provides an arrangement for essentially 
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performing, in the case of a triangular mesh, operations that generally provide the inverse of Loop's 
surface subdivision methodology, and, in the case of a quadrilateral mesh, operations that generally 
provide the inverse of the Catmull-Clark surface subdivision methodology. 

Loop's surface subdivision methodology will be described in connection with FIG. 2. 
Generally, in Loop's surface subdivision methodology, each triangular face in the original mesh is 
split into a plurality of sub face, the sub faces defining the finer level mesh. The vertices of the finer 
level mesh are positioned using weighted averages of the positions of the vertices in the original 
mesh. More specifically, and with reference to FIG. 2, a mesh 20 at a fine "j" level is depicted that 
includes a vertex v q 21(0) located at position d(q) 9 and a plurality "K" of surrounding points v q (l) 
through v q (K) (generally, v q (k)) 21(1) through 21(K) (generally identified by reference numeral 
21(k)) located at positions at respective positions d(l) through c^K), the points v q (k) 21(k) being 
those points in the mesh 20 that are connected to vertex v q by edges. ("K," the number of points 
connected to vertex v q will sometimes be referred to as the vertex's "valence. 1 ') Denoting the set of 
the indices of the points v q (k) in the mesh that are connected to vertex v q by N(q j), a mesh at the 
next higher fineness level "j+1" corresponds to vertex v' q 22(0) located at position c i+1 (q) and a set 
of surrounding points v f q (l) 22(1) at located at positions d +l Q) and connected to the vertex v' q 22(0). 
The finer level mesh is constructed by providing the vertex v' q 22(0) at a position d +1 (q) that is 
determined by 




(1) 



and points 



v' q (l) 22(1) at positions that are determined by 



c J+l U) = \[c j (q)+c j {k)]+\[c J {k-\)+c j {k+\)],k=\,...,K (2), 



where the weighting factor a(K) is given by 
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f (In 
3+ 2 cos — 
v K ) 

8 



(3). 



1 

2 Equations (1) and (2) are applied by selecting each point in the mesh at level "j" as a vertex, to 

3 provide points for the mesh at level "j+1 ." Loop's surface subdivision methodology can be applied 

4 recursively to provide a mesh at any desired degree of fineness. It will be appreciated that equation 
f!5 (2) will provide that the point v' q (l) 22(1) that is created when the methodology is applied to point 
;:|> v. 21(0) as the vertex will be in the same position when the methodology is applied to point v (k) 

21 (k) as the vertex. 

p| It will be appreciated that, the Loop surface subdivision methodology provides 

m 

■xr • 

3 9 (i) one new point, at location c i+1 (l) in the finer ("j+1") level mesh, for, and generally 

P 

ftf) somewhere between, each pair of points in the coarse ("j") level mesh, and 

i\ (ii) a relocated point, at location c i+I (q) in the finer ("j+1") level mesh for, and generally 

§L2 somewhat near, each point in the coarse ("j") level mesh, 

13 which are interconnected to form the finer ("j+1) level mesh. Accordingly, in performing the 

14 inverse, that is, processing a finer ("j+1") level mesh to generate a coarse ("j") level mesh, the 

15 arrangement processes the points in the finer ("j+1") level mesh to remove, for each pair of 

1 6 connected points in the finer ("j+1 ") level mesh, and determines a position for the remaining point, 

1 7 which will be present in the coarse ("j ") level mesh, based on the positions of surrounding points in 

18 the finer ("j+1") level mesh. The operations performed by the arrangement in determining the 

19 position of the point that is to remain in the coarse ("j") level mesh will be described below. 

20 The Catmull-Clark subdivision methodology produces smooth surfaces using a small number 

2 1 of neighboring vertices. The Catmull-Clark surface subdivision methodology will be described in 
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1 connection with FIGS. 3A through 3E. With reference to FIG. 3A, that FIG. depicts a mesh 30 

2 consisting of four quadrilaterals 31(1) through 3 1(4). Each quadrilateral is referred to as a face. The 

3 quadrilaterals are defined by points 32(0), which is common to all of the quadrilaterals 31(1) through 

4 31(4), and other points 32(1) through 32(8). Generally, the Catmull-Clark surface subdivision 

5 methodology is performed in a series of iterations, including 

6 (i) a face point generation iteration, in which, for each face, a face point 33(1) through 33(4) 

7 is generated as the average of the points 32(p) defining the respective face (reference FIG. 3B); 

8 (ii) an edge point generation iteration in which, for each edge, a new edge point 34(1) 

9 through 34(4) is generated as the average of the midpoint of the original edge with the two new face 
10 points of the faces adjacent to the edge (reference FIG. 3C; the face points generated during the face 
M point generation iteration are shown as circles); 

2 (iii) a vertex point iteration in which a vertex point 35 is generated in relation to the positions 

f i 3 

ip of the new face points generated in iteration (i) for faces that are adj acent to the original vertex point, 
the positions of the midpoints of the original edges that are incident on the original vertex point, the 

45 position of the original vertex point, and the valence of the original vertex point (reference FIG. 3D; 

S$ the face and edge points generated during the face point and edge point generation iterations are 
shown as circles); and 

18 (iv) a mesh connection step (reference FIG. 3E) in which: 

1 9 (a) each new face point is connected to the edge points of the edges defining the original 

20 face; and 

21 (b) each new vertex point is connected to the new edge points of all original edges 

22 incident on the original vertex point. 

23 More specifically, in the Catmull-Clark subdivision methodology, denoting, for a vertex v q 

24 at location c*(q), 

25 (i) the set of indices of neighboring points d(l) that are connected thereto (for example, for 

26 point 32(0), point 32(2) and 32(8)) by N e (q,j), and 
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(ii) the set of indices of second-order neighboring points that are opposite to vertex d(q) with 
respect to the level "j" faces (for example, for point 32(0), point 32(1)) that are incident with vertex 
ctobyN&j) 

a mesh at the next higher fineness level "j+1 " is constructed as follows. In the face point generation 
iteration, the face points are generated and located at positions c i+1 (m i ) determined as follows 



m t e N f (g,j+ l),/„/, +1 € N e (gj),n e N f (qj) 



In the edge point generation iteration, the edge points are generated and located at positions d (\ { ) 
determined as follows: 

^'M = iH^+^+^^+^w), (s) 

h e N e (q,j+ l),*, e W e (4,7), ro, e N f (q,j+ l) 



In the vertex point generation iteration, the new vertex points are generated and located at positions 
as follows: 



l^NXj+lqlm^N^j+lq) 



In terms of only points c*, that is substituting for the last term in equation (6), 



C J 



V 4KJ 2 K leNe[qj) 4K meN/[qj) 



(7). 
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1 

2 Generally, for quadrilateral meshes, the arrangement makes use of the Catmull-Clark 

3 methodology, except for the case of K=3 in equation (6). In that case, the arrangement makes use 

4 of 

c J+{ (q) = {l-Sr)c J (q)+4r^ !>(/)+ 4^ £ c J+ \m) (8), 

A lzN e {q,j) A m&N f (q,j+\) 



6 where y = — . Equation (8) would correspond to the Catmull-Clark methodology (equation (6)) 
38 

Si 3 1 

1:J with Y = , except that, for K=3, r = — instead of — . In terms of only points d (compare 

w 1 4K 1 38 12 

LI! 

ru 

equation (7)) 

L c*+\q) = {l-7r)c J (q)+6r^ T,c J {l) + r ^ Z c j {m) (9 ). 

y i 

u.i 

10 As is apparent from the above discussion, the Catmull-Clark surface subdivision 

1 1 methodology applied to a quadrilateral mesh at the coarse ("j") level will produce a quadrilateral 

1 2 mesh at the finer ("j+1") level having the same number of points as the mesh at the coarse ("j ") level, 

13 but the points in the finer ("j+1") level may be at different positions than in the coarse ("j") level. 

14 Accordingly, in performing the inverse, that is, processing a finer ( n j+l ") level mesh to generate a 

15 coarse ("j") level mesh, the arrangement processes the points in the finer ("j+1") level mesh 

16 determine a position for each point in the finer ("j+1") level mesh based on the positions of 

17 surrounding points in the finer ("j+1") level mesh. The operations performed by the arrangement 

1 8 in determining the position of the point in the coarse ("j") level mesh will be described below. 

19 As noted above, the invention provides an arrangement for generating a coarse ("j) level 

20 mesh from a finer ("j+1") level mesh. In generating a coarse mesh from a fine mesh, an indicator 
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function is used locally in connection with each point that is being processed as a vertex to determine 
whether, around the respective point, the fine mesh can be derived by subdivision of a coarser mesh 
using a subdivision inverse filter. If so, the value of the indicator function, which will be described 
below, will be "zero" or close to "zero." On the other hand, if the fine mesh cannot be derived by 
subdivision of a coarser mesh in the region around the respective vertex using a subdivision inverse 
filter, the value of the indicator function will not be close to "zero," and in that case the coarse mesh 
will be generated using an approximate solution of a least-squares problem. 

The indicator function is constructed for each point, selected as a vertex v f q , using the value 
of the Laplacian L(k,j+1) of the positions of the point and of the other points in its neighborhood in 
the finer (j+1) l eve l mesh. For triangular meshes, in the one-dimensional case, for example, for a 
vertex v' q on a boundary, crease line or the like in the finer ("j+1 ") level mesh, the Laplacian for the 
vertex is 



L(kJ+ l) = \[c J+ '(k- 1) + c J+l (k+ 1)]- c J+l (k) 



(10), 



where c j+1 (k) is the position of the vertex v' q in the finer ("j+1) level mesh for which the Laplacian 
is being generated, and c i+1 (k-l) and c j+1 (k+l) are the positions of the neighboring points v q (k-l ) and 
v' q (k+l) in the fine (j+1) level mesh. In the two-dimensional case, that is, in the case of a vertex v ? q 
other than on a boundary, crease line or the like, for a regular vertex, that is, for a vertex for which 
the valence "K" is equal to "six," the Laplacian is 

\( ) 

L( k >J +l ) = 7 2> +, (/) -c^{k) (11), 



where the factor "1/6" reflects the valence of the vertex. Also in the two-dimensional case, for an 
irregular vertex, that is, for a vertex for which the valence "K" is not equal to "six," which may be 
the case at, for example, edges of the mesh, the Laplacian is 
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L(kJ+l) = p 



K 



(12), 



1 

2 where 

3 3+Sa(K) 

P = - 3(-5 + Sa(K)) (13) ' 

I'M where a(K) is defined in equation (3) above. 

j^f) For quadrilateral meshes, in the one-dimensional case, that is, for a vertex v' q on a boundary, 

crease line, or the like in the finer ("j+1") level mesh, the Laplacian L(k,j+1) for the vertex v' q is 

UJ 

E_ 5 

f L(k,j+\) = ^(c J+] (k-\)+c J+i (k+l))-c J+l (k) (14) 

Ul 

where c i+1 (k) is the position of the vertex v' for which the Laplacian is being generated in the finer 

Mi H 

("j+1 ") level mesh, and (^(k-l) and c i+I (k+l) are the positions of the neighboring points in the finer 

10 (j+1) level mesh. In the two-dimensional case, that is, in the case of a vertex v' q other than on a 

1 1 boundary, crease line or the like, two Laplacians are used, namely a Laplacian L e (k,j+1) using the 

12 position d +l (k) of the vertex V and the positions of the points whose indices are in the set N e (q,j), 

13 and a Laplacian L^kj+l) using the position d +1 (k) and the positions of the vertices whose indices 

14 are in the set Nf(q,j), as follows 



\l€N e (kJ+\) 



- C j+ \k) (15) 



16 and 
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A ( h ... A 



- C 



7 + 1 



(16). 



As will be described below, the arrangement makes use of values of both Laplacians L e (k,j+1) and 
Lf(kj+1), and so both values will be generated for each vertex. 

As noted above, the arrangement makes use of an indicator function. For triangular meshes, 
the indicator function i(q) for a point as vertex v' q is given by, in the one-dimensional case 

i{q) = L(q) - (L{q - l) + l{q + l)) + \ (l{q - l) + l{q + 2)) (17) 
and, for the two-dimensional case 

i(q) = -6L{q)+2 %L{!)-2 X L(l) + ^ L(l) (18)> 

leN(\,0)(q) leN(2,0)(q) /e^(2,l)( 9 ) 

where the indices for the sums in equation (18) are in the neighborhood of the vertex for which the 
indicator function is being evaluated are as depicted in FIG. 4. 

In generating the values of the indicator function for quadrilateral meshes, the arrangement 
makes use of a linear combination L ind (kj+1) of the Laplacians L e (kj+1) and Lf(kj+1) described 
above (equations (15) and (16), respectively), namely 

L ind {kJ+\) = 2^L t {kJ+\)+^L f {kJ+l) (19) 

for vertices for which valence "K" is not equal to "three." For vertices for which the valence "K" 
does equal "three," L ind (kj+1) takes the form 
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*w(*../+l) = -yl i (tj+l)+21 / (tj+l) (20). 

2 For regular vertices, that is, for vertices for which the valence "K" equals "four/' equation (19) 

3 becomes 

4 LZ = -2L,(kJ+l)+L f (kJ+l) (21). 

5 With this background, the indicator function i(k) for quadrilateral meshes is given by (omitting the 

6 index "j+1," since the indicator function will always be evaluated for vertices in the fine (j+1) level 
Ql mesh) 

-~v 

ti\ 

I *(*) = fu*)- EUO^ 2XM + 

:r ( 22 )> 

= Z /eJV(3,l)(*) Z /etf(3,2)(Jt) H /etf(4,2)(Jfc) 

where the vertex neighborhoods for the respective sums are as shown in FIG. 5. 

3 

* i 0 As noted above, the fine-to-coarse mesh generating arrangement generates a coarse ("j ") level 

I I mesh from a finer ("j+1") level mesh. For a triangular mesh, the fine-to-coarse mesh generating 

12 arrangement determines, for each vertex v' q located at position c i+1 (k) in the finer ("j+1") level mesh 

13 for which a vertex v q is to appear in the coarse ("j") level mesh, a position d(k) for the vertex v q in 

14 accordance with 

15 c J {k) = c j+ \k) + X L(kJ+ l) (23), 

16 where X is a parameter whose value are determined as described below and L(k,j+1) is the value of 

17 the Laplacian generated for the vertex v' q in the finer (j+1) level mesh. As noted above, not all of 
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the vertices in the finer ("j+l ") level mesh will have corresponding vertices in the coarse ("j") level 



for a selected small e, then for equation (23), the position d(k) of the corresponding vertex v q in the 
coarse ("j") level mesh can be taken as the position d + \k) of the vertex v ! q . In that case, neither the 
parameter A,, nor the indicator function i(q) (equations (17) and (18)), need to be generated. 

If the value of the Laplacian L(k,j+1) is not small, that is, if the condition in equation (24) 
does not hold, the manner in which parameter X for use in equation (23) is to be determined will 
depend on the value of the indicator function, generated as described above in connection with 
equations (17) and (18) and whether the vertex is on a boundary, crease line or the like, a regular 
vertex (that is, a vertex v' whose valence "K" is "six") or an irregular vertex (that is, a vertex v' q 
whose valence "K" is other than "six' 1 ). If the value of the indicator function is zero, or close to zero, 
for the vertex v' q , a subdivision-inverse filter is used in which 

(i) if the vertex is on a boundary, crease line or the like 



mesh. If the value of the Laplacian L(k j+l ) generated for a vertex 
is relatively small, that is, if 



v' q in the finer ("j+l") level mesh 




(24) 



X = -1 



(25), 




(iii) if the vertex v' q is irregular, that is, if its valence "K" is not equal to "six," and is not on 
a boundary, crease line or the like 
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X = 



8a(K) 
-5+ %a{K) 



(27), 



where a(K) is as defined above in connection with equation (3). 

On the other hand, if the value of the indicator function i(q) is not zero or close to zero, a 
least-squares optimization is used in which 

(i) if the vertex is on a boundary, crease line or the like 



X = 



L{k) 



b\ D L{k) + \b\ D [L{k- l) + L(k+ l)) 



(28), 



ID 12 W 23 

where b 0 = - — and b x = - — . 



(ii) if the vertex v' q is regular, that is, if its valence "K" is equal to "six," and the vertex is not 
on a boundary, crease line or the like 



X = 



1 

L(k) 



b r <* L(k) + \b[ es ££(/) 



(29), 



61 „ 14403 1 

where K eg = and b, = ; it will be appreciated that the factor " — " in equation 

0 5720 1 5720 6 

(29) reflects the value of the valence "K," namely, "six." 

(iii) if the vertex v ? q is irregular, that is, if its valence "K" is not equal to "six," and the vertex 
is not on the boundary, crease line or the like 
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X = 



L(k) 



leN(kJ+\) 



(30), 



where 



2 (5- Sa(K)) (l 4647 K- 391848 a{K) + 391848 a{kf) 
0 " 715(3+ Sa(k))(256+ 41K- 512 a(k) + 256 a{k) 2 ) 



and 



16(-553L^- 24521 a( K) + 24521 a(K) 2 ) 
1 715(256+ 41^- 512 a(K)+ 25 a{K) 2 ) 



where, for both equations (31) and (32), "K" corresponds to the valence of the irregular vertex, and 
a(K) is as defined in equation (3). 

For a quadrilateral mesh, for each vertex v' q located at position c i+1 (k) in the finer (j+1) level 
mesh, a vertex v q is provided in the coarser (j) level mesh at a position c^k) that is determined by 

c J (k) = c J ^ l (k)+A l LXkJ+l)+A 2 L f (k 9 j+l) (33), 



where X } and are parameters whose values are determined as described below, and L e (k,j+1) and 
L^kj+l) are determined as described above in connection with equations (14) through (16). In the 
one-dimensional case, that is, for vertices on a boundary, crease line or the like, since only Laplacian 
L e is used, only one parameter X x need be generated for use in equation (33). On the other hand for 
vertices that are not on a boundary, crease line, or the like, since both Laplacians L e and L f are used, 
both parameters X, and ^ are generated for use in equation (33). As with the triangular mesh, if the 
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values of L e (kj+1) and L^kJ+l) are relatively small, that is, if for a vertex v' q in the fine ("j+1") 
level mesh 



max (|l e (*»•/ + l)|, ^(£,7 + l)|) < e 



(34) 



for a selected small e, then for equation (34), the position c'(k) of vertex v q in the coarse ("j") level 
mesh can be taken as the position c j+l (k) of vertex v' q . In that case, neither the parameters X,, and ^ 
nor the indicator function i(kj+l) (equation (22)), need to be generated. 

If the condition in equation (34) does not hold, the values of parameter X, and ^ for use in 
equation (33) will depend on the value of the indicator function. For a vertex v' q that falls on a 
boundary, crease line or the like, the value of the parameter X l will be generated in the same manner 
as in a triangular mesh (reference equations (25) and (28) above). On the other hand, for a vertex 
v' q that does not fall on a boundary, crease line or the like, if the value of the indicator function 
i(kj+l) generated therefor is near zero, that is, if 



i{k,j+l)\ < £-max^L e (k,j+l)\,\L f (k,j+l)\j 



(35) 



a subdivision-inverse filter is used in which, if the valence "K" of the vertex not equal to "three," 



X 
X 



1 K- 3 

1 

2 " K-3 



(36), 



and, if the valence "K" for the vertex is equal to "three," 

X j = - 8, X 2 = 2 (37). 
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On the other hand, if, for a vertex v' q , the value of the indicator function i(k j+1 ) is not near zero, that 
is, if the condition in inequality (35) is not satisfied, 



1 



6£4(*,y+l)+iaff I £.(/,./+ l) 

A leN(kJ+\) 



(38), 



and 



L f {k,j+l) 



^ leN f {k,j+\) 



(39), 



which reflect a least squares optimization. If vertex v' q is regular, that is, if the vertex's valence "K" 
is "four," 



ICC 

&10 



ICC 



ucc 
°20 



rcc 
°2\ 



9946871 
4862025 
1024 



405 

1644032 

972405 
1338874 

972405 



(40). 



On the other hand, if the vertex v' is irregular, that is, if the vertex's valence "K" is other than "four," 
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rCC 

b \o 



ICC 



ICC 

b 20 



ICC 

b 2\ = 



162307143936- 92746939392 if - 8924282387 K 3 
4862025 (12544 - 14336 K + 4096 K 2 + 901 K 3 ) 

1024(2793728- 1596416/:- 244001 K 3 ) 
99225(12544- 14336 K + 4096 A: 2 + 901 K 3 ) 

512 (-113305472+ 64745984^+ 17391149 K 3 ) 
4862025(12544- 14336^+ 4096 K 2 + 901 K 3 ) 

4(8660934688- 4949105536 1876158821 K 3 ) 
4862025 (12544 - 14336 K + 4096 K 2 + 901 K 3 ) 



(41). 



The following are approximations for equations (40) that are used in one embodiment of the 
arrangement: 



7 CC 

°\o = 


- 2.04583 


ICC 

b n = 


- 2.5284 


ICC 

b 20 = 


1.69069 


ICC 

b 2\ = 


- 1.37687 



(42), 



and the following are approximations for equations (41) that are used in the same embodiment: 
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1623.07 - 927.469 K - 89.2428 K 3 
609.892 - 697.02 K + 199.149 K 2 + 43.8068 K 3 

28.6078 - 16.3473 K - 2.49857 K 3 
12.4468- 14.2249*:+ 4.06426 K 2 + 0.894017 K 3 

(43). 

-580.124+ 33 1.499 # + 89.0427 K 3 
609.892 - 697.02 K + 199.149 K 2 + 43.8068 K 3 

346.437 - 197.964 K - 75.0464 K 3 
609.892 - 697.02 K + 199.149 K 2 + 43.8068 K 3 

1 

Q. With this background, the operations performed by the fine-to-coarse mesh generating 

rip arrangement in generating a coarse ("j") level mesh from a finer ("j+1 ") level mesh will be described 

1% in connection with the flow charts in FIGS. 6 and 7. The flow chart in FIG. 6 depicts the operations 

U5 performed in connection with a triangular mesh and the flow chart in FIG. 7 depicts the operations 

(Pp performed in connection with a quadrangular mesh. With reference initially to FIG. 6, the 

*_7 arrangement will first perform a series of steps to generate the Laplacian (equations (10) through 

LIS (12)) for each of the points in the triangular finer ("j+1") level mesh. Initially, the arrangement will 

3* ""2 

select a point in the finer (j+1) level mesh as vertex (step 100) and determine whether the point is 

IE? 

Bb regular, that is, whether its valence "K" is six, and does not fall on a boundary, crease line or the like 

1 1 (step 101). If the arrangement makes a positive determination in step 101, it proceeds to step 102 

12 in which it generates the Laplacian for the point selected in step 100 in accordance with equation 

13 (11). 

14 On the other hand, if the arrangement makes a negative determination in step 101, it will 

1 5 determine whether the point is irregular, that is, whether its valence "K" is less than six and does not 

16 fall on a boundary, crease line or the like (step 103). If the arrangement makes a positive 

17 determination in step 103, it proceeds to step 104 in which it generates the Laplacian for the point 

18 selected in step 100 in accordance with equation (12). 



010 



ICC 



ICC 

b 20 



ICC 

b 2\ 
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1 Finally, if the arrangement makes a negative determination in step 103, which will be the 

2 case if the point falls on a boundary, crease line or the like, it proceeds to step 105 in which it 

3 generates the Laplacian for the point selected in step 100 in accordance with equation (10). 

4 Following step 102, 104 or 105, the arrangement will determine whether all of the points 

5 have been selected in step 1 00 (step 1 06). If the arrangement makes a negative determination in step 

6 106, it will not have generated Laplacians for all of the points in the finer (' j+1 ") level mesh, and, 

7 in that case, it will return to step 100 to select another point and perform respective steps 101 

8 through 105 in connection therewith. The arrangement performs step 100 and respective steps 101 

9 through 1 05 through a plurality of iterations, in each iteration selecting and generating the Laplacian 
10 for each of the points in the finer ("j+1") level mesh. 

= s 
si: 

it After all of the points have been selected and Laplacians generated therefor, the arrangement 

[ii will make a positive determination in step 106, at which point it will sequence to a series of steps 

;tj$ in which the arrangement determines the positions c'(q) for points in the finer ("j+1") level mesh that 
are to be retained in the coarse ("j") level mesh. In those operations, the arrangement will select a 

Hi 

J 5 point in the finer ("j+1 ") level mesh (step 1 1 0) and determine whether the point is to be retained in 

JtfS the coarse ("j") level mesh (step 111). If the arrangement makes a negative determination in step 

£117 1 1 1 , the point is not to be retained in the coarse ("j") level mesh. In that case, the arrangement will 

Jtfe sequence to a step 1 12 in which it determines whether it has selected all of the points in the finer 

rV9 ( M j+1 ") level mesh. If the arrangement makes a negative determination in step 1 12, it will return to 

20 step 1 10 to select another point in the finer ("j+1") level mesh. On the other hand, if the arrangement 

2 1 makes a positive determination in step 1 1 2, which will be the case if it has processed all of the points 

22 in the finer ("j+1") level mesh, it will exit (step 113). The arrangement can make use of any of a 

23 number of methodologies in determining whether a point in the finer ("j+1") level mesh will be 

24 retained in the coarse ("j") level mesh. For example, given a suitable point indexing methodology, 

25 as will be apparent to those skilled in the art, the arrangement can determine whether a point in the 

26 finer ("j+1") level mesh will be retained in the coarse ("j") level mesh. 

27 Returning to step 1 1 1 , if the arrangement makes a positive determination in that step, which 

28 will be the case if the selected point in the finer ("j+1 ") level mesh will be retained, the arrangement 

29 sequences to step 115, in which it determines whether the value of the Laplacian for the point 
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selected in step 1 10 is less than "e" (reference inequality (24)). If the arrangement makes a positive 
determination in step 115, it determines the position c"(q) of the point in the coarse ("j") level mesh 
as the position c i+1 (q) of the corresponding point in the finer ("j+1") level mesh (step 116). 
Following step 1 16, the arrangement will sequence to step 1 12 to determine whether all of the points 
in the finer ("j-H") level mesh have been processed, and, if not, return to step 1 10 to select another 
point, and, if so, exit. 

Returning to step 1 1 5, if the arrangement makes a negative determination in that step, it will 
sequence to step 1 17, in which it will determine whether the point is regular, that is, whether its 
valence "K" is six, and does not fall on a boundary, crease line or the like. If the arrangement makes 
a positive determination in that step, it will generate a value for the indicator function (reference 
equation (18)) (step 118). The arrangement will then determine whether the value of the indicator 
function is zero or close to zero (step 119). If the arrangement makes a positive determination in 
step 1 19, it will provide a value for the parameter X as described above in connection with equation 
(26) (step 120). On the other hand, if the arrangement makes a negative determination in step 119, 
it will provide a value for the parameter X as described above in connection with equation (29), 
provided that the value of the parameter X does not fall below the value specified in equation (26) 
(step 121). 

Returning to step 1 17, if the arrangement makes a negative determination in that step, it will 
sequence to step 122 to determine whether the point is irregular, that is, whether its valence "K" is 
not equal to six, and the point does not fall on a boundary, crease line or the like. If the arrangement 
makes a positive determination in that step, it will generate a value for the indicator function 
(reference equation (18)) (step 123). The arrangement will then determine whether the value of the 
indicator function is zero or close to zero (step 124). If the arrangement makes a positive 
determination in step 124, it will provide a value for the parameter X as described above in 
connection with equation (27) (step 125). On the other hand, if the arrangement makes a negative 
determination in step 124, it will provide a value for the parameter X as described above in 
connection with equation (30), provided that the value of the parameter X does not fall below the 
value specified in equation (26) (step 126). 
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1 Returning to step 1 22, if the arrangement makes a negative determination in that step, it will 

2 sequence to step 127 to determine whether the point falls on a boundary, crease line or the like. If 

3 the arrangement makes a positive determination in that step, it will generate a value for the indicator 

4 function (reference equation (17)) (step 128). The arrangement will then determine whether the 

5 value of the indicator function is zero or close to zero (step 129). If the arrangement makes a 

6 positive determination in step 129, it will provide a value for the parameter X as described above in 

7 connection with equation (25) (step 130). On the other hand, if the arrangement makes a negative 

8 determination in step 129, it will provide a value for the parameter X as described above in 

9 connection with equation (28), provided that the value of the parameter X does not fall below the 
1 0 value specified in equation (25) (step 131). 

U Following step 120, 121, 125, 126, 130, or 131 , the arrangement will sequence to step 132, 

fi| in which it determines whether the value of the parameter X generated in the respective step was 

13 greater than "one," and, if so, limit the maximum value of parameter X to "one" (step 133). 

£3- Thereafter, the arrangement uses equation (23) to determine the position c i+1 (q) of the point in the 
coarse ("j") level mesh (step 134). 

S 

€h Following step 134, the arrangement the arrangement will sequence to step 1 1 2 to determine 

ftf whether all of the points in the finer ("j+1") level mesh have been processed, and, if not, return to 

±p step 1 10 to select another point, and, if so, exit. 

1 9 The flow chart in FIG. 7 depicts the operations performed in connection with a quadrangular 

20 mesh. With reference to FIG. 7, the arrangement will first perform a series of steps to generate the 

2 1 Laplacians (equations (14) through (16)) for each of the points in the quadrilateral finer ("j+1 ") level 

22 mesh. Initially, the arrangement will select a point in the finer (j+1 ) level mesh as vertex (step 1 50) 

23 and determine whether the point does not fall on a boundary, crease line or the like (step 151). If 

24 the arrangement makes a positive determination in step 151, it proceeds to step 152 in which it 

25 generates the Laplacian for the point selected in step 1 50 in accordance with equations (15) and ( 1 6). 

26 On the other hand, if the arrangement makes a negative determination in step 151, which will 

27 be the case if the point falls on a boundary, crease line or the like, it proceeds to step 1 53 in which 

28 it generates the Laplacian for the point selected in step 150 in accordance with equation (14). 
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1 Following step 152 or 153, the arrangement will determine whether all of the points have 

2 been selected in step 1 50 (step 1 54). If the arrangement makes a negative determination in step 1 54, 

3 it will not have generated Laplacians for all of the points in the finer ("j+1 ") level mesh, and, in that 

4 case, it will return to step 1 50 to select another point and perform respective steps 1 5 1 through 153 

5 in connection therewith. The arrangement performs step 1 50 and respective steps 151 through 1 53 

6 through a plurality of iterations, in each iteration selecting and generating the Laplacian for each of 

7 the points in the finer ("j+1") level mesh. 

8 After all of the points have been selected and Laplacians generated therefor, the arrangement 

9 will make a positive determination in step 154, at which point it will sequence to a series of steps 
1 0 in which the arrangement determines the positions d(q) for respective points in the coarse ("j ") level 
lA mesh using the positions c" +1 (q) of the points in the finer ("j+1") level mesh and the values of 
|2 parameters to be generated and the Laplacians generated in steps 150 through 153. Initially, the 
p arrangement will select a point in the finer ("j+1 ") level mesh (step 1 60) and determine whether the 

values of both Laplacians L e (k,j+1) and L^kj+l) generated for the point selected in step 160 are less 

S'f than "e" (reference inequality (34)) (step 161). If the arrangement makes a positive determination 

*L6 in step 16 1 , it determines the position d(q) of the point in the coarse ("j ") level mesh as the position 

UF ^ +l (q) of the same point in the finer ("j+1") level mesh (step 162). Following step 162, the 

H$ arrangement will sequence to step 163 to determine whether all of the points in the finer ("j+1") level 

Wt) mesh have been processed, and, if not, return to step 1 60 to select another point. On the other hand, 

20 if the arrangement makes a positive determination in step 163, it will exit (step 164). 

2 1 Returning to step 1 6 1 , if the arrangement makes a negative determination in that step, it will 

22 sequence to step 1 65, in which it determines whether the determine whether the point is regular, that 

23 is, whether its valence "K" is four, and does not fall on a boundary, crease line or the like. If the 

24 arrangement makes a positive determination in step 165, it will generate a value for the indicator 

25 function (reference equation (22)) (step 166) and determine whether the value of the indicator 

26 function is near zero (reference inequality (35)) (step 167). If the arrangement makes a positive 

27 determination in step 1 67, it will generate values for A,, and as described above in connection with 

28 equations (36) (step 168). On the other hand, if the arrangement makes a negative determination in 

29 step 167, it will generate values for X, and as described above in connection with equations (38)- 
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1 (40), in the process ensuring that the values of X l and 7^ are within predetermined limits (-4< ^ <; 1/2, 

2 1/4*^*1) (step 169). 

3 Returning to step 1 65, if the arrangement makes a negative determination in that step, it will 

4 sequence to step 170, in which it determines whether the point is irregular, that is, whether its 

5 valence "K" is other than four, and does not fall on a boundary, crease line or the like. If the 

6 arrangement makes a positive determination in step 1 70, it will generate a value for the indicator 

7 function (reference equation (22)) (step 171) and determine whether the value of the indicator 

8 function is near zero (reference inequality (35)) (step 172). If the arrangement makes a positive 

9 determination in step 172, it determine whether the valence "K" of the point is "three" (step 173). 
1 0 If the arrangement makes a positive determination in step 173, it will generate values for X } and ^ 
KI as described above in connection with equation (37) (step 174). On the other hand, if the 

arrangement makes a negative determination in step 173, it will generate values for X x and X^ as 
described above in connection with equation (36) (step 175). Returning to step 172, if the 

4% arrangement makes a negative determination in step 1 72, which will be the case if the value for the 

Hi 

55 indicator function as generated in step 171 is not near zero, it will generate values for X } and as 

!L6 described above in connection with equations (38)-(40), in the process ensuring that the values of 

^ ( 4 1 1 1 ^ 

£IJ7 X, and X l are within predetermined limits y - — — - < X x <— , — < A 2 < — — -J (step 176). 



1 8 Returning to step 1 70, if the arrangement makes a negative determination in that step, it will 

19 sequence to step 177 to determine whether the point falls on a boundary, crease line or the like. If 

20 the arrangement makes a positive determination in that step, it will generate a value for the indicator 

21 function (reference equation (17)) (step 178). The arrangement will then determine whether the 

22 value of the indicator function is zero or near zero (reference inequality (35) above) (step 179). If 

23 the arrangement makes a positive determination in step 1 79, it will provide a value for the parameter 

24 X as described above in connection with equation (25) (step 180). On the other hand, if the 

25 arrangement makes a negative determination in step 179, it will provide a value for the parameter 

26 X as described above in connection with equation (28) (step 181). 
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1 Following step 168, 169, 174, 175, 180, or 181, the arrangement will sequence to step 182, 

2 in which it uses equation (33) and values for parameters \ and generated in the previous steps to 

3 determine the position c* +1 (q) of the point in the coarse ("j") level mesh. Following step 182, the 

4 arrangement will sequence to step 1 63 to determine whether all of the points in the finer ("j+ 1 ") level 

5 mesh have been processed, and, if not, return to step 1 60 to select another point, and, of so, sequence 

6 to step 1 64 to exit. 

7 The arrangement provides a number of advantages. In particular, the invention provides an 

8 arrangement for generating a coarse ("j ") level mesh representing a surface, from a finer ("j + 1" ) level 

9 mesh surface representation. The arrangement makes use of an indicator function to provide an 
1 0 indication as to which of several methodologies is be used at respective ones of the points in the finer 
13 Cj + 1") level mesh, including a subdivision-inverse filter methodology and a least-squares 
|(2 optimization methodology. 

3 ?? 

4h It will be appreciated that numerous changes and modifications may be made to the 

SCft arrangement as described herein. For example, it will be appreciated that, although the arrangement 

sat. 

T5 has been described in connection with meshes of triangular and quadrangular faces, it will be 

ff 5 appreciated that the arrangement may find utility in connection with meshes having faces of different 

Ui 

fW polygonal structures. In addition, although equations for generating the indicator functions and 
parameters may differ than the specific ones described herein. 

iL-i- 

1 9 It will be appreciated that a system in accordance with the invention can be constructed in 

20 whole or in part from special purpose hardware or a general purpose computer system, or any 

2 1 combination thereof, any portion of which may be controlled by a suitable program. Any program 

22 may in whole or in part comprise part of or be stored on the system in a conventional manner, or it 

23 may in whole or in part be provided in to the system over a network or other mechanism for 

24 transferring information in a conventional manner. In addition, it will be appreciated that the system 

25 may be operated and/or otherwise controlled by means of information provided by an operator using 

26 operator input elements (not shown) which may be connected directly to the system or which may 

27 transfer the information to the system over a network or other mechanism for transferring 

28 information in a conventional manner. 
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1 The foregoing description has been limited to a specific embodiment of this invention. It will 

2 be apparent, however, that various variations and modifications may be made to the invention, with 

3 the attainment of some or all of the advantages of the invention. It is the object of the appended 

4 claims to cover these and such other variations and modifications as come within the true spirit and 

5 scope of the invention. 

6 What is claimed as new and desired to be secured by Letters Patent of the United States is: 
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